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Abstract 

We report new measurements of the decays B + — > p + ^f, B° — > p°7 and B° — > W7 us- 
ing a data sample of 657 x 10 6 B meson pairs accumulated with the Belle detector at the 
KEKB e + e- collider. We measure branching fractions B(B+ p + j) = (8.7^1?;?) x HT 7 , 
B(B° -»• p° 7 ) = (7.8 ±i;J ±?;g) x 1(T 7 and B(B° -»• w 7 ) = (4.0 t 1 ^ ± 1.3) x HT 7 . We also report the 
isospin asymmetry A(/ry) = — 0.48lg'^g ig'og; and the first measurement of the direct CP-violating 
asymmetry Acp(B + — > p + j) = — 0.11±0.32±0.09, where the first and second errors are statistical 
and systematic, respectively. 

PACS numbers: 11.30.Hv, 13.40.Hq, 14.65.Fy, 14.40.Nd 
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The b — > d-y process, which proceeds via a loop diagram (Fig.[TJa)) in the Standard Model 
(SM), provides a valuable tool to search for physics beyond the SM, since the loop diagram 
may also involve virtual heavy non-SM particles ll). The process has been observed in the 
exclusive modes B — >■ fry and B — > coy by Belle jlj and Babar [jjj]. Branching fractions for 
these modes have been used to constrain the ratio of Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements Q |Vtd/T4 a |; a non-SM effect may be observed as a deviation of |V^/V ts | 
from the expectation based on measurements of other CKM matrix elements and unitarity 
of the matrix An additional contribution from an annihilation diagram (Fig. EJb)) may 
induce a direct CP-violating asymmetry in B + — > p + 7, and an isospin asymmetry between 
B — > fry modes; the latter can be used to constrain the CKM unitarity triangle angle 03 j6[. 
These quantities are also sensitive to physics beyond the SM Q . In this paper, we report new 
measurements of the B — > fry and B — > wy processes using a data sample of (657 ± 9) x 10 6 
B meson pairs accumulated at the T(4S) resonance. With a data sample almost twice as 
large and an improved analysis procedure, these results supersede those in (2)]. 

(a) loop diagram 
w 

FIG. 1: (a) Loop diagram for 6—7-^7 and (b) annihilation diagram, which contributes only to 
B + -»■ p+7. 

The data are obtained in e + e~ annihilation at the KEKB energy-asymmetric (3.5 on 
8 GeV) collider jsj and collected with the Belle detector (9)]. The Belle detector in- 
cludes a silicon vertex detector (SVD), a central drift chamber (CDC), aerogel threshold 
Cherenkov counters (ACC), time-of-flight (TOF) scintillation counters, and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals located inside a 1.5 T superconducting 
solenoid coil. An iron flux-return located outside of the coil is instrumented to identify K\ 
and muons (KLM). 

We reconstruct three signal modes, B + — > p + 7, B° — > p°7 and B° — > ouj, and two control 
samples, B + — > K* + y and B° — > K* 0/ y. Charge-conjugate modes are implicitly included 
unless otherwise stated. The following decay modes are used to reconstruct the intermediate 
states: p + ->■ 7r + 7r°, p° ->■ 7r + 7r _ , lu -)■ 7r + 7r _ 7r , K* + ->■ K + n°, K*° ->■ K + n~, and 7r° ->■ 77. 

Photon candidates are reconstructed from ECL energy clusters having a photon-like shape 
and no associated charged track. A photon with an T(4S") center-of-mass (cm.) energy (E*) 
in the range [1.8, 3.4] GeV is selected as the primary photon candidate. Photons detected by 
the endcap ECL, which were excluded in the previous analysis, are also used. To suppress 
backgrounds from tt /t] — > 77 decays, we apply a veto algorithm based on the likelihood 
calculated for every photon pair consisting of the primary photon and another photon. We 
also reject the primary photon candidate if the ratio of the energy in the 3x3 crystal array, 
centered on the crystal with the maximum energy, to that in the 5x5 array is less than 
0.95. 

Neutral pions are formed from photon pairs with invariant masses within ±16 MeV/c 2 
(~3<r) of the 7T° mass. We require the energy of each photon to be greater than 50 MeV, 
and the cosine of the angle between the two photons in the laboratory frame to be greater 
than 0.58 (0.40) for the 7r° from p + (u). The photon momenta are then recalculated with a 
7T° mass constraint. 



(b) annihilation diagram 
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Charged pions and kaons are selected from tracks in the CDC and SVD. Each track is 
required to have a transverse momentum greater than 100 MeV/c and a distance of closest 
approach to the interaction point within 0.5 cm in radius and within ±3.0 cm along the 
positron beam (z-) axis. We use a likelihood ratio C w /(jC w + £#) < 0.3 for pions and > 0.6 
for kaons, where the pion and kaon likelihoods C n and Ck are determined from ACC, TOF 
and CDC information. The criteria have efficiencies of 86%, 87% and 89% for a pion from 
p + , p° and co, respectively; the misidentification probability for a kaon is 8.3% (8.5%) for 
p + Kaons for K* candidates are selected with an efficiency of 85%. Invariant masses 

for the p, lo and K* candidates are required to be within windows of [640,890] MeV/c 2 , 
[760,800] MeV/c 2 , and [820,970] MeV/c 2 , respectively. 

Candidate B mesons are reconstructed by combining a p or oj candidate with the pri- 
mary photon and calculating two variables: the beam-energy constrained mass M^ c = 



(i?b cam /c 2 ) 2 — \pg/c\ 2 , and the energy difference AE = E* B — E£ eam . Here, p£ and E* B are 
the cm. momentum and energy of the B candidate, and E^ eam is the cm. beam energy. To 
improve resolution, the magnitude of the photon momentum is replaced by (E^ eam — E*^)/c 
when the momentum p B is calculated. 

To optimize the event selection, we study Monte Carlo (MC) events in a signal box defined 
as 5.273 GeV/c 2 < M bc < 5.285 GeV/c 2 a nd \AE\ < 0.1 GeV. For each signal mode, we 
choose selection criteria to maximize Ns/y/Ns + Nb, where N$ and Nb are the expected 
signal and the sum of the background yields. 

The dominant background arises from continuum events (e + e~ — > qqij), q = u,d,s,c), 
where a random combination of a p or uj candidate with a photon forms a B candidate. 
We suppress this using a Fisher discriminant (J 7 ) calculated from modified Fox- Wolfram 



moments [10] and other variables, i.e., the cosine of the polar angle (cos6 B ) of the B direc- 
tion, the distance along the z-axis (Az) between the signal vertex and that of the rest of 
the event and, in addition, M^+^-^o and Dalitz plot variables for the wy mode. For each of 
these quantities, we construct likelihood distributions for signal and continuum events. The 
distributions are determined from MC samples. 

From these likelihood distributions we form likelihoods L s and C c for signal and contin- 
uum background, respectively. In addition, we use a flavor-tagging quality variable r that 
indicates the level of confidence in the 5-flavor determination as described in Ref. In 
the (r, 7Z) plane defined by the tagging quality r and the likelihood ratio 1Z = C s / (C a + C c ), 
signal tends to populate the edges at r = 1 and 7Z = 1, while continuum preferentially 
populates the edges at r = and 7Z = 0. We divide the events into six bins of r (two 
bins between and 0.5, and four bins between 0.5 and 1) and determine the minimum 7Z 
requirement for each bin. In the p +r y mode, if the tagging-side flavor is the same as that of 
the signal side, we assign the events to the lowest bin < r < 0.25. The TZ criteria reject 
98% of continuum background while retaining 35%, 51% and 43% of the p°7, p +r y and u-y 
signals, respectively. For the K* +/ y (K* 0r y) mode, we use the criteria for the p + j (p°j) mode. 

We consider the following backgrounds from B decays: B — > K*j, other B — > 
processes, decays with a 7r°/r) (B — > pre , utc°, prj and oorj), other charmless hadronic B 
decays, and b — > c decay modes. The B — > K*^ background can mimic the B — > p7 signal if 
the kaon from the K* is misidentified as a pion. To suppress B — > K*^ events we calculate 
Mk-k, where the kaon mass is assigned to the charged pion candidate; for p°7, the lower of 
the two Mk-k values is taken (misassignment tends to give a higher Mk-k)- For the p + 7 mode, 
we reject the candidate if Mkk < 0.90 GeV/c 2 , while for the p°7 mode we use Mkk in the 
fit procedure to extract the signal (note: Mkk > 0.92 GeV/c 2 is required when optimizing 
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selection criteria). The b — > 37 modes (B — > K*^ and other B — > X s ^y decays) contribute 
to the background when the p and u candidates are formed from random combinations of 
particles. Decays with a 7r°/?7 can mimic the signal if one of its daughter photons is not 
detected. To suppress this background, we reject the candidate if | cos^hdl > 0.80, 0.75 and 
0.83 for the p + 7, p°7 and wy mode, respectively, where the helicity angle ^hei is the angle 
between the tt + track (the normal to the u decay plane) and the B momentum vector in 
the p ibj) rest frame. We study large MC samples and find no other distinctive hadronic B 
decay background sources. 

The reconstruction efficiency for each mode is defined as the fraction of the signal remain- 
ing after all selection criteria are applied, where the signal yield is determined from a fit to 
the sum of the signal and continuum MC samples using the procedure described below. We 
take the pion identification efficiency from a data sample of D* + — > D°tt + , D° — > K~n + . 
The total efficiencies are listed in Table HI The systematic error in the efficiency is the 
quadratic sum of the following contributions, estimated using control samples: the uncer- 
tainty in the photon detection efficiency (2.4%) as measured in radiative Bhabha events; 
charged tracking efficiency (1.0% per track) from partially reconstructed D* + — > D°7T + , 
D° — > Kg7r + TT~ , Kg — > 7r + (7r~); charged pion and kaon identification (0.5 to 0.6% per track) 
from D* + — > D°7i + , D° — > K~7t + ; neutral pion detection (4.6%) from 77 decays to 77, 
7T + 7r~7r° and 37r°; and TZ-r and 7i°/r] veto requirements (2.0 to 8.4%) from B — > Dtt~ with 
D° ->• K-n + , K° s u and D+ K-n + n + . 

We perform an unbinned extended maximum likelihood fit to Mb c and AE (and M^k 
for the p°7 mode) for candidates satisfying |Ai?| < 0.5 GeV and Mb c > 5.2 GeV/c 2 . The 
fit is performed individually for the three p'y/wy signal modes and the two K*^ modes. We 
describe the events in the fit region using the sum of probability density functions for the 
signal, continuum, K*^ (for the p7 modes only), and other background hypotheses. We 
use the distributions of MC events in histograms to model the M^ C -AE shapes of B decay 
background components and the Mk-k shapes for all components. 

The signal distribution for the p°7 and K* ^ modes is modeled as the product of a Crystal 
Ball lineshape [l2| in AE to reproduce the asymmetric ECL energy response, a Gaussian in 
Mbc, and an MC histogram distribution for Mx n . For the p + 7, K* +/ y and modes, we use 
the product of a Crystal Ball lineshape for AE and another Crystal Ball lineshape for M^ c . 
The signal parameters of M^ c and AE shapes for K*^ modes are determined from fitting 
the data; for the P7/W7 modes, they are taken from MC and calibrated using the data/MC 
difference of the fits to the K* +r j and K* 0/ y samples for the modes with and without a neutral 
pion, respectively. 

The continuum background component is modeled as the product of a linear function 



in AE, an ARGUS function 13] in Mb c , and, for p°7, an MC histogram for Mk-k- The 
continuum shape parameters and normalizations are mode dependent and allowed to float. 

There is significant K* 0/ y background in the p°7 sample. This background is modeled by 
the product of a two-dimensional M^ C -AE histogram and an Mk-k histogram. Similarly, the 
K* +r y background for p + 7 is modeled by a two-dimensional M^ C -AE histogram. In both 
cases, the AE peak position is shifted from the p7 signal peak; this offset is determined from 
fitting the MC histogram shape to a K*^ data sample in which the pion mass is assigned 
to kaons. The same K*^ sample together with the known kaon to pion misidentification 
probability is also used to determine the size of the K*^ background. 

Other B decays are considered as an additional background component when we extract 
the signal yield. The levels of these backgrounds are fixed using known branching fractions 
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FIG. 2: Projections of the fit results to M bc (in |A#| < 0.1 GeV and 0.92 GeV/c 2 < M K7T ), AE (in 
5.273 GeV/c 2 < M bc < 5.285 GeV/c 2 and 0.92 GeV/c 2 < M K7r ), and for 5° -> p°y, M K?r . Curves 
show the signal (dashed, red), continuum (dot-dot-dashed, blue), B — > K*y (dotted, magenta), 
other backgrounds (dash-dotted, green), and the total fit result (solid). 



or upper limits 14 



The systematic error in the signal yield due to the fitting procedure is estimated by 
varying each of the fixed parameters by ±lcr and then taking the quadratic sum of the 
deviations in the branching fraction from the nominal value. The varied parameters are the 
signal shape parameters, branching fractions of the background components, AE shift of 
the K*j component, and the kaon to pion misidentification probability determined from a 
control sample. The results of the fits are shown in Fig. [2] and listed in Table [I] 

The systematic error in the branching fraction has contributions from the efficiency, fit- 
ting, and the number of B meson pairs; we add these together in quadrature. The significance 



is defined as y — 2 ln(£ /^max), where £ max (£o) is the value of the likelihood function when 
the signal yield is floated (set to zero). To include systematic uncertainty, the likelihood 
function from the fit is convolved with a Gaussian systematic error function. 

Tabled also lists combined branching fractions, which are calculated from the products of 
likelihoods from individual fits. We combine p + 7 and p°7 modes (referred to as fry) and three 
P7 and wy modes (referred to as (p, w) 7 ) assuming a single branching fraction B(B — > fry) = 
B(B (p,u/)7) = B{B + p+ 7 ) = 2 x ^B(B° ->• p° 7 ) = 2 x ^B(B° -»• uj) 0, Efj, 

I L ^ ^ I — I 

where = 1.071±0.009 [17|. The results are consistent with the previous measurements [2|, 
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TABLE I: Yield, significance with systematic uncertainty, efficiency, and branching fraction (B) 
for each mode. The first and second errors in the yield and B are statistical and systematic, 
respectively. The sub-decay B(u — > tt + tt^tt°) is included for the ^7 mode. 

Mode Yield Significance Efficiency (%) B (10~ 7 ) 

B+ -> p+ 7 45.8 tfd tH 3~3 8.03 ±0.59 8.7lf?i°i 

B o^ p o 7 757 +16.8 +5.1 50 14.81 ±0.95 7.8^1^ 

B°^uj~f 17.5 tfjt^ 2.6 6.58 ±0.76 4.0 ±^±1.3 

B-^pj 5.8 12.1 ±1.2 

B^(p,u;)j 6.2 11.4 ± 2.0 



0] and have smaller errors. They are also in agreement with SM predictions p, [la, 18 

The ratios of the branching fractions of the B -+ pj/oo^ modes to those of the B — >■ 7^*7 



modes can be related to |Vfd/Vt s | p, lla]. We calculate the ratios from likelihood curves of 
individual fits to the B — > pj/u>j and B — > K* r y samples. Systematic errors that do not 
cancel in the ratio are convolved into the likelihoods. We find, 



B(B° -+ p° 7 ) 
B(B° -+ K* "/) 
B(B -+ p 7 ) 
B(S -+ ^7) 
B(S^(p,a;)7) 



n nonfi +0.0045 +0.0014 /-. \ 

U.UZUU _o.0043 -0.0016; \ L ) 

n n^no +0.0060 +0.0026 

U.UOUZ _ .0055 -0.0028' K^l 

0.0284 ± 0.0050 (3) 



B[B -> K*7) 

where the first and second errors are statistical and systematic, respectively. 



Using the prescription in Ref. |6|, Eq. |3]for example gives \Vu/V ts \ = 0.195 taoi9( ex P-) ^ 
0.015(th.). This is consistent with determinations from B® mixing [14], which involve box 
diagrams rather than penguin loops. We also find B(B + — > K* +/ ~f) = (384 ± 17) x 10~ 7 
and B(B° — > K* 0, ~f) = (378 ± 8) x 10~ 7 (statistical error only), in agreement with the world 
average. 

From Table HI we calculate the isospin asymmetry A(/ry) = ^—^-B(B + — > p +r y)/B(B° — > 
p°7) — 1 and find 

A(p 7 ) = -o.48±8S±8:SI- (4) 

The result is in agreement with the previous measurement jjjj], and is only marginally con- 



lne result is m agreement witn trie pr( 
sistent with the SM expectations [g 1 , 0] • 

We also calculate the direct CP-violating asymmetry A C p(B + — > p + j) = [iV(p~7) — 
iV(p + 7)]/[iV(p _ 7) + iV(p + 7)] using a simultaneous fit to B + — > p +r y and B~ — > p~7 data 
samples. We consider systematic errors due to the fitting procedure, asymmetries in the 
backgrounds, and possible detector bias estimated using aB + Dn control sample. We use 
the measured asymmetries (TiJ for B + — > K* +/ ~f, p + vr°, p + r] and B — > X s 7 and assume up 
to 100% asymmetry for other charmless hadronic B decays. We find 

A CP (B + p+7) = -0.11 ±0.32 ±0.09. (5) 

The result is consistent with the SM predictions p 1 , [l6[ . 
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In conclusion, we present a new measurement of branching fractions for B — > p7 and 
B — > ury, a measurement of the isospin asymmetry, and the first measurement of the direct 
CP-violating asymmetry for B + — > p + 7. The results are consistent with SM predictions. 
We improve the experimental precision on |V^/VL| determined from penguin loops, finding 
good agreement with the value determined from box diagrams 14 . 
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